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Abstract 

The inverse spectral transform for the Zakharov-Shabat equation on 
the semi-line x > is reconsidered as a Hilbert problem. The boundary 
data induce an essential singularity as k —* oo to one of the basic solutions. 
Then solving the inverse problem means solving a Hilbert problem with 
particular prescribed behavior. It is demonstrated that the direct and 
inverse problems are solved in a consistent way as soon as the spectral 
transform vanishes as 0(1/ k) at infinity in the whole upper half plane 
(where it may possess single poles) and is continuous and bounded on the 
real fc-axis. The method is applied to stimulated Raman scattering and 
sine-Gordon (light cone) for which it is demonstrated that time evolution 
conserves the properties of the spectral transform. 

1 Introduction 

Since the discovery of the inverse spectral transform (1ST) in 1967 |[] to solve 
nonlinear evolutions (Korteveg-de Vries and later nonlinear Schrodinger j2|, then 
sine-Gordon, etc... [|3|) with Cauchy conditions on the infinite space (and van- 
ishing asymptotic boundary values), the application of the method to boundary 
value problems on semi-infinite (or finite) intervals has been the object of in- 
tense research with few general results, which is reminiscent from the fact that 
even for linear partial differential equations, boundary value problem much dif- 
fer from one another. It would be useless to quote all papers that has dealt with 
boundary value problems for nonlinear evolutions, only the most recent ones 
will be cited here, see Jj|-[p5| where more complete reference will be found. 

Taking for instance the nonlinear Schrodinger equation (NLS) for a field 
q(x,t), saying that it is solvable by 1ST on the infinite line actually means 
that, given Dirichlet condition q(x,0) on the line x £ (— oo, +oo) and Cauchy 
conditions g(±oo, t) — and q x (±oo, t) — for any t > 0, the method 1ST allows 
to construct the solution q(x, t) at any t for all x. Although Dirichlet conditions 
at x = ±00, namely g(±oo, t) — 0, constitute with q(x, 0) a well-posed problem, 
the solution of NLS with 1ST requires the adjunction of the Neuman condition 
q x (±oo,t) = 0. This simply results from the fact that the Lax pair contains 



1 



cxplicitely the derivative q x . The point is that this supplementary constraint is 
compatible with the rest of the boundary values and remains compatible at any 
time. This stability is a property of the spectral transform itself that conserves 
the class of function (potentials) chosen at t — 0. For instance if q(x, 0) is in the 
Schwartz space (the function and all its derivatives are continuous and vanish 
faster than any power of x as x — > oo), then the solution q(x, t) obtained by 1ST 
is also in the Schwartz space which makes the method work. 

Such is not anymore the case on finite (or semi-infinite) interval. Indeed, 
while the constraint q x (±oo,t) — is not restrictive when we already have 
g(±oo,<) = 0, it is not so on a finite domain where imposing the data of both 
the function and its derivative produce an over-constrained problem. Then, al- 
though NLS keeps many of its integrability properties on the finite interval (with 
mixed Dirichlet-Neumann conditions on the semi-axis) (J] [0, and Backlund 
transformation can be used to determine explicit solvable subclasses |J, the 
general problem is still unsolved. 

A detailed study of the Dirichlet boundary problem (data of q(x,0), q(0,t) 
and q(oo, t) = 0) has been performed in |l0| where it has been shown that 
the problem of the determination of q x (Q, t) from g(0, i) is transformed into the 
problem of determining some missing spectral datum. This approach allows also 
to construct classes of solutions but not to solve the problem in general. 

Another point of view is presented in where the method is based on 
the scattering problem on a line in the (x, i)-plane in the Gelfand-Levitan- 
Marchenko formulation. This is applied to the Korteveg-de Vries equation and 
also provides much insight in the boundary- value problem but shows that indeed 
the problem is still open. 

An interesting approach has been proposed in ll^ j where both operators of 
the Lax pair are treated as joint spectral problems. There the boundary-value 
problem becomes a Dirichlet condition on a closed domain, namely data of 
q{x,0), q(0,t), q(oo,t) = and q(t, oo) = 0. Applied to NLS (defocusing case), 
this method provides the solution as a system of coupled integral equations 
which is then difficult to use. But this approach allows for particular solutions, 



see also 12 , and for large time asymptotic behaviors. Moreover it has been 
succesfuly applied to sine-Gordon in the light cone (SG) (llj for which the Lax 
pair is compatible with the boundary values, and then extended to boundary 
values on arbitrary convex domain 113] . The method has also proved being very 
useful to treat linear partial differential equations with moving boundaries Jl5| . 
The method furnishes the solution under the form of a set of Cauchy-Green 
integral equations. 

There is another case where the general solution has been obtained by 1ST: 
the stimulated Raman scattering equations (SRS) on the semi-axis |l6) (the 
sharp-line version on the finite interval has been treated in p7| by the method 
of fll3[). The boundary- value problem there is an open-end problem and the 
stability of direct/inverse spectral transformation allows indeed to get the gen- 
eral solution on the semi-axis (as well as the finite interval case) which has been 
succesfully checked on numerical simulations Jlifl . The method works also for 
the discrete integrable version of SRS 119] and it is based on the approach of 
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the dispersion relation |2fl] where the principal Lax operator is taken as the 
spectral operator and the auxiliary Lax operator is the tool to define altogether 
the dispersion relation and the time-evolution of the spectral transform. The 
advantage of such a method is that, in the physically relevant case of a medium 
initially at rest, it allows to get all the information on the nonlinear Fourier spec- 
trum explicitely from the input radiation fields. An interesting consequence is 
the possibility to build literaly the spectral transform out of measurements and 
compare it to the theory, which has been done successfuly on numerical simu- 
lations EsJ of the finite interval, as well as experimental data plj with infinite 
line damped model. The method has also been extended to treat perturbation 
of SRS [p2[ , confirming the validity of approximations used for comparison with 
experiments in pl[ . 

As a general result, it has been shown that the dispersion relation will depend 
explicitely not only on the boundary values (as for infinite line |p3| ) but also on 
the spectral transform itself hence producing nonlinear time evolution of spectral 
data (of Riccati type) (Hi. Interesting enough, such an evolution induces the 
motion of the poles (discrete eigenvalues of the spectral problem, associated to 
solitons) and thus results in creation (and anihilation) of solitons due to input 
boundaries. The same property has been found for the karpman-Kaup system 
[ p4| with a remarkable adequation to numerical simulations ]25| . 

The spectral problem underlying stimulated Raman scattering and sine- 
Gordon is the Zakharov-Shabat system on the semi-axis with prescribed po- 
tentials values in a; = 0. We show here that the solution can be written out 
of a system integral equations of Cauchy type where the contour of integration 
differs from the one working in the infinite line case. This results from the 
presence of an essential sigularity in the asymptotic behaviors of one of the Jost 
solutions, a problem that has been underestimated in the litterature so far. 

Then we shall begin with solving the direct and inverse spectral problems 
in sections || and ||. In doing this we shall discover that the reconstructed 
solutions (Jost solutions) must obey a sequence of properties, named PI to P9, 
that are central to the method. This is then applied to SRS in section |] and 
to sine-Gordon in section |^. In both cases it is demonstrated that the time 
evolution of the spectral transform is such that it conserves the constraints 
necessary to ensure the 9 fundamental properties of the Jost solutions. In order 
to illustrate the method for the sine-Gordon, we consider the explicitely solvable 
case of a piecewise constant boundary datum which generically stands for the 
discretization of an arbitrary boundary value. 



2 The direct problem 
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2.1 Basic solutions 

The time t being an external parameter it is omitted in the whole section. We 
consider the Zakharov-Shabat spectral problem [|| on the semi-line 

x>0 : 4> x =ik{<P,cx 3 } + Q{x)4> , Q=(° r J) 
with given boundary values 

<l( x )\ x =0+ = lo , r(x)\ x=0+ = r . (2.2) 
Jost solutions |2(| can be defined as the solution of the integral equations 

( i>ti(k,x) \ _ ( 1 \ / ffdy q(y)^ti(k,v) \ (9 ^ 

1, ^(k,x) ) {0 ) + {j Q x dy r(y)i>+(k,y)e 2ik (*-ri ) ' W 

( <t>UKx) \_ ( \ / -£°dy q {y)4>Uk,y)e-^-^ \ 

\ 4>Uk,x) ) \ 1 ) + \ f x dyr(y)cj)t 2 (k,y) J ' [ZA) 

\ _ ( 1 \ / f*dy q(y)(j> 21 (k,y) \ (e) ,x 
fc(k,x) J~{0 ) + \~jrdy r{v)<t>n{k,v)eW-v) ) ^ 

^ 2 (k,x) \_f 0\ f J Q x dy q(y)^ 2 (k,y)e- 2ik ^ \ 
fe(fc,x) J \l J + \ fidy r(y^- 2 (k,y) ) ' ^ 

2.2 Boundary behaviors 

The so-called physical solution 

*(k,x) = , (2.7) 



obeys the bounds 



<f(MUoo=( V r+ ) ' ^ 

where the reflection coefficients p (k) are defined by reading the above integral 
equations at x = + , namely 



P + (k) = - dy q(y)<f>l 2 (k,y)e 2ik y 
Jo 

p-(k) = - dy r{y)^ 1 (k,y)e- 2lk y , (2.10) 



(i 
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and the transmission coefficients (fc) by making the limit as x — » oo i.e. 



T+(fc) = 1 + 

r"(fc) = 1 + 



dy r(y)<pf 2 (k,y) , 

) 



(2.11) 



By direct computation with (2.1) we demonstrate that the determinant of any 
solution is independent on x and hence can be calculated at x = or at x — > oo. 
Applied to <& it gives 



1 — p p = T T 



which is called the unitarity relation. 
The rest of the solution, namely 



V(k,x) = (iP+,^) 



(2.12) 



(2.13) 



is used to solve the inverse problem and we call it the intermediate solution. It 
obeys the following bounds for k £ M 



i o 
o 1 



l/r H 



..,2'ikx 



p /t- 



1/t- 



(2.14) 
(2.15) 



(the boundary behavior as x — > oo is actually obtained by using the relations 
that follow). The intermediate solution is related to the physical solution by 
the following relations for k G R 



4>t - V^ 2 



,.,2ikx 



P 4) 2 



_ -2ifcc „+„/,+ 



(2.16) 



which can be proved in different ways, for instance by the computation of the 
differences out of (2.4) fl2.6p and the comparison the resulting integral equation 
with the one for the r.h.s. of ( 2.16| ). A simpler but equivalent method is to 
express the vector solution on the basis ipi and e 2lkx Tp2 ■ The coefficients of 
this expansion are then calculated by comparing the boundary values in x = + . 
Another proof goes with expressing the solution ^> in terms of $ as 



1 p+(k) 
p-{k) 1 



a ik<r 3 x 



(2.17) 



(computed by matching the values in x = + ) which in particular readily pro- 
vides the behavior (2.15) from (2.9). 
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2.3 Analytical properties 

Using standard method |2(J one obtains that, for potentials vanishing fast 
enough as x — ► oo, the vector ip^ is analytic in the upper half fc-plane and 
ip 2 i n the lower half. Next, the vector cf> 2 is meromorphic function of k in the 
upper half-plane with N + simple poles £;+ which are the zeroes of the analytic 
function l/r + (fc). Indeed from (2.15) we have 



#21 



(2.18) 



with tp 21 analytic. Moreover the vector 



can be written from (2.4) as the solution of 



<pj 2 (k,x) 
^ 2 (A;,x) 



/* dy q(y)4>i 2 (k iy )e-™(*-rt 
Jo dy r{y)4>i 2 (k,y) 



(2.19) 



(2.20) 



which shows that (j) 2 is analytic in k and thus that <p 2 has poles where the 
analytic function l/r + (fc) has zeroes. 

Similarly, 07/ is meromorphic in the lower half-plane with N~ poles k, 
zeroes of l/r~(fc). Consequently, from their very definitions |TTo|) (fTO]) , t 
and p^ have meromorphic extensions in their respective half-planes ±Im(fc) > 
where they possess the simple poles k„ (the bound states locations '). In 
particular we have the relation, obtained simply by taking the residues on (2.17) 



n ' 
± 



Res{0+} = p+V 1 +(fc+)exp[-2i^ 



= Res{p+} 

k+ 



R-es{</>i } = P n i } 2 (K ) exp[2i/c„ x] , p n = Rcs{p } 

k^ k~ 

The spectral data is then constituted by the set 

S = {Im(fc) > , p + {k) ; Im(fc) < , p~ (fc)} , 



(2.21) 
(2.22) 

(2.23) 



which is complete if we prove that it allows to reconstruct the potentials r(x) 
and q(x) with q(0 + ) = qo and r(0 + ) = ro. 



2.4 Large k behaviors 

In order to solve the Hilbert problem ( 2.16 ) with the singularities k^ in the 
complex plane, one must calculate first the behaviour of the solutions <&(£;, x) 
and ^ (k, x) as |fc — > oo. Note that we follow Mushkhclishvili |27 to refer to 
( [2.16 ) as a Hilbert problem. 
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By using repeatidly the integration by part we obtain from the integral 



equations ( |2.4| )(2.6) the following regular asymptotic behavior of the physical 
solution 



1 



) 



+ — ( ~ ft rq 3 

2ik \ -r J Q rq 

- Jo M - r %Jo r «] ~ qx + q Jo r l 
(2ik) 2 V ~r x + r J£ rq - J* [rq y - rq J£ rq] 



O 



fc 3 



(2.24) 



which has to be understood in Im(k) < for the first vector (f> 1 and in Im(fc) > 
for 4%. 

The same procedure applied to the intermediate solution shows that it pos- 
sess an essential singularity at k — oo on the real axis as indeed we get 



{^1,^2 
1 



1 



L + —(~ft rq £ 

2ik \ -r J Q rq ) ' 2ik \ r e 

- In l r v1 - r %ft r l} -<lx + qJ X rq 
(2ik) 2 V ~r x + r J* rq - J* [rq v - rq J Q V rq] 

1 ( (qe 2zkx - q )r (q f* rq + q' )e~ 2tkx 
{2ik) 2 \ (r Jo rq + r' Q )e 2tkx [ re - 2ikx - r )q 





2iki 



-2ikx 



(2.25) 



Consequently, when qo and ro do not vanish, the intermediate solution has not 
the same behavior as the physical solution, in contrast with the infinite line 
case. In particular the formulation of the solution of the Hilbert problem will 
differ from the one of the infinite line. 

Note for future use that from the definitions ( 2.10 ) of p^- follow the large k 
expansions 



p+{k) = ik qo ~Wk) 



2 l'o 



1 



p-(k) = 



1 



1 



i 



2ik 

where we have defined 



ro 



(2ik) 



2 '0 



{2ikf W-rlqo) + G^) 



% = 



dq{x) 



dx 



„ _ d 2 q{x) 
% ~ dx 2 



(2.26) 



(2.27) 



x=0+ 



so as for r(x). 



3 Solution of the Hilbert problem 

3.1 Statement of the problem 

Let be given the set of spectral data 

S = {Im(fc) > , p + (k) ; Im(fc) < , p~{k)} 



(3.1) 
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where p + (k) (resp. p~{k)) is continuous and bounded on the real axis and 
meromorphic in the upper (resp. lower) half plane with a hnite number of 
single poles fc+ (resp. k~) and related residues p+ (resp. p~). Moreover they 
obey in their respective half-planes 

p+{k) = <k qo+ ° { ¥ h p ~ {k) = -'k r0+o{ T* ) (3 - 2) 

where qo and tq are given numbers. The problem is to construct the two matrices 
&(k, x) and ^(k, x) solution of the Hilbert problem (2.16) obeying the following 
set of properties. 

P. 1 : Alternative expression of the Hilbert problem 

$ = *M , M = e~ lkcJ3X ( ^_ P ^ e lka3X , x > . 

P.2 : The vector (resp. tp^) i s analytic in Im(fc) > (resp. Im(fc) < 0). 
The vector </>J (resp. is meromorphic in Im(fc) > (resp. Im(fc) < 0) where 
it possess N + (resp. N~) simple poles fc+ (resp. k~) with residues 

Res{<^} = p+ipf(k+)exp[-2ik+x] , Res{<^} = Pn^iK ) exp[2i/c~ x] . 

fe+ 

P.3 : The functions ^(k,x) and $>(k,x) obey the bounds 

•cm>-(j ;). «(*.o)-(/ w « + f> 

P.4 : the function fy(k,x) — 1 vanishes for x < 0. 

P.5 : the function $(fc, x) has a Laurent expansion with a pole of zero order 
at k = oo: 

$(fc, x) = 1 + \<S> {1) {x) + 0(l/k 2 ) 
k 

(in the lower half-plane for the first vector, in the upper half-plane for the 
second). 

P.6 : the function ^(fc, x) has a Laurent expansion with an essential singu- 
larity at k = oo on the real axis 

* = 1 + y^^(x) + l*W e 2ifc^3X + ... 

k k 

P.7 : There exists an off-diagonal matrix Q(x) such that for x > 

V x = ik[V,<T 3 ] + Q(x)$ , 

which is called the Zakharov-Shabat spectral equation. 
P. 8 : The function Q obeys 

Q(z)L= 0+ = ( r ° o o ) ' G(*)L=o- = . Q( x ) ■ 

P. 9 : The above defined function Q(x) is also given by 

Q(x)=i[a 3 ,$W(x)} . 
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3.2 The solution 

We prove now that the solution of the Hilbert problem ( 2.1 6| ) obeying the prop- 
erties P.1-P.9 above is given by the following equations. First $(fc, x) is given 
by the explicit expressions 

/ 1 \ 1 r+oc+iO i-k 

< M " ( J ) " 3EF L+„ —k ^(VNrft.) 

"EfI^^-^" 1 ' lm(fc)<0 (3.3) 

n=l n 

«<*•*> -(!)+S5 JCT ATI 



^+(fc+,x)e- 2ifc " :E , lm(fc)>0 (3.4) 



n=l 



where the function ^(k,x) solves the following Cauchy-Green integral system 
of equations 

/ 1 \ i p+oo—iO j\ 



Pn „;,- , ;.- ..,,2/7... 
* ^ "/ 



$^(*£,a:)e**» a , Im(fc) > (3.5) 

/■,' . -- -■ 

n=l 



n \ 1 /•+oo+i0 ji 

U » 1 ' " A ^ ,-2^,,,+ , 



N+ + 



iif(ki,x)e- 2lk " x , lm(fc)<0 (3.6) 



In contrast with what occurs in the one dimensional case x € (— oo, +oo), here 
the integrals do not run on the real axis but on contours along the real axis in 
the upper or lower half-plane as indicated. We shall see that this apparently 
slight difference has a lot of consequences and allows in particular to reconstruct 
the very solutions $ and ^ defined in the direct problem. 

It is worth stressing here that we do not adress the problem of the complete 
characterization of the spectral data p ± (fc), namely to find the necessary and 
sufficient conditions which would ensure a unique solution the the above Cauchy- 
Green intergral equations (such a problem has been worked out in jl^]). However 
we shall work with a set of sufficient conditions on p ± (fc) ensuring that the 
reconstructed potential does coincide with the one we started with. In particular 
for what concerns solution of nonlinear evolutions, we will demonstrate that a 
set of actual spectral data (constructed from a given potential at time 0) evolve 
in time such as to conserve this set of sufficient conditions. 

We proceed now step by step to demonstrate that the expressions (3.3)-(3.6) 
do produce functions <& and \& that verify properties P.l to P.9. 
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Property P.l : For real values of k the integral equations ( |3.5[ ) (3.6) can be 
written by virtue of Cauchy theorem (by assumption p ± (fc) have no poles on 
the real axis) 



1 



+oo+i0 



dX 

iO A — 



p-(X)e 2iXx ^(X,x) 



N~ 



(k)e 2 * kx ^(k,x)-Y, 



Pn 



n=l kn 



-ip 2 (k n , x)e 



2ik„ x 



k G 



(3.7) 



\ 1 f+°°- i0 dX 



^ {K x) =[l)+^ J__ iQ — k ^ (A)e — V J" (A, x) 

N+ + 

- p+(k) e - 2ik ^+(k, x)-J2 7-T^T^(k+, x)e~ 2ik ^ , keR. (3.8) 

n=l kn k 

Hence follow the Hilbert relations 

4>t - ^ = ~e 2ikx p-^2 . 4>t ~ 4>2 = e- 2lkx P + ^f , (3.9) 
and by rearrangement of the vectors into matrices, P.l is proved. 

Property P.2 : The analytical properties of the various vectors follow from 
their very expressions. Indeed, as the set of single poles of $ is explicitely 
displayed, one has only to prove that the Cauchy-like integrals produce holo- 
morphic functions |27j . Such is the case if the integrant is a continuous bounded 
function of A el vanishing as A — > oo. This actualy results from the assumed 
spectral data p ± (X) and from the large k e M properties of ty(k,x) deduced 
from the integral equations (|3.5| ) and ( |3.6| ) as the property P. 6 that follows. 



Property P.3 : The equations ( |3.5|) ( p.q ) evaluated at x — {x is a parame- 
ter) lead to the system 



N~ 

-E^fc-O), (3-10) 
n=l Kn K 

N+ + 
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Now, as we assume a unique solution to the integral equations, it is sufficient to 
verify that ^(k, 0) = 1 is indeed solution of the above system. For ^{k, 0) we 
have then the equation 

^-(i)-i£rA^(J)-|i^(!)- 

which by contour integration in the lower half plane, and thanks to the properties 
of p~ , easily gives ip^(k,0) = (1,0) T (remember that k lives in the upper half- 
plane). The same procedure applies to the second vector and gives ^{k, 0) = 
(0,1) T 

Next, by using exactly the same method with $(fc, 0), and taking care of the 
pole A = k, we do obtain the property P.3 (note that one could also use P.l to 
derive the boundary value of one of the functions from the other one). 

Pro pert y P.4 : For x < 0, thanks to the factor e 2lXx appearing in the expres- 
sion ( |3.5| ), we can perform the integration in the lower half-plane for if)^. The 
contribution of the poles exactly cancels out and ip^(k,x) = (1,0) T for x < 0. 
The same procedure holds for ip% which proves P.4. 



Property P. 5 : T he e xpansion of Q(k,x) at large k readily results from the 
epressions ( |3.3| ) and (3.4) where the Cauchy integral behave well. Indeed, taking 
e.g. (f>i , we integrate above the real axis where, for x > 0, the exponential is 
bounded. 

We obtain in particular for the matrix $W(a;) = (<j>i , $2) 



.'"--L/ p -(\)e 2 ^^(\,x) + Y,Pn4>2(K,x)e 2lKx , (3-12) 



1 2ilT 



— 00+iO 



1 /•+(»— tO N 

'" ' p + (A) e - 2 ^+(A,x) + ^p+^+(fc r t,a;)e- 2 ^ . (3.13) 



2 2iir 



°°- i0 n=l 



Property P.6 : The expansion of ^(k,x) at large k S M. cannot be ob- 



tained directly from the integral equations (3.5) and (3.6) (one sees on (3.5) 
that exp[2iAa:] is integrated for x > on a contour lying in the lower half A- 
plane), but rather form the versions ( |3.7| ) and (3.S) where the Cauchy integrals 
behave well. Then the same procedure as for P.5 applies to get the expansion 
given in P.6 (it is actually easier to use P.5 and P.l together with the expan- 
sion of the spectral data p ± {k) to demonstrate P.6). We obtain in particular, 
according to the notations in P.6, 



V 2i 



-3a 

' 2i 



(3.14) 



which matches the expressions (2.25) (and can be pursued to next orders). 
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Property P.7 : This is a fundamental property in the process of solving 
the inverse problem as indeed it gives the expression of the potential from the 
solution of the Cauchy-Green integral equations ( 3.£ ) (3.6), in other words in 
terms of the spectral data p ± (k). 

To perform the proof of P.7 it is first more convenient to rewrite the integral 
equations as follows 

^ = (o)-il^" e2i ^ 2 "' (3 - 15) 
o \ j_ r 

2in J c X - k 



K = ( i + 2* L x^-""^ ■ (3 - 16) 



where C+ (resp. C_) is a smooth contour in the upper (resp. lower) half-plane 
extending from — oo to oo and passing over all poles k„ (resp. under all poles 
fc~). Note that tp^ is written for Im(fc) > and hence there is no contribution 
of the pole A = k in the lower half-plane. The same holds reversely for tp^ . In 
short these integral equations read 

1 f dX 

*(fe, x) = 1 + — / — -¥(A, x)R(X, x) , (3.17) 



where 

R(k,x) = e~ ik ^ ( _ p ° {k) P+ Q (fc) ) e ik °°* , (3.18) 

and where the contour of integration has to be understood as being C_ for the 
first vector of the product ^R, and C + for the second. 

Thanks to the particular dependence in x of R(k,x) it is easy to compute 
\& x + ik[a 3 , ty] and by algebraic manipulations to get 

V x + ik[<r a , *] - - J d\[a 3 , *R] 
1 f d\ 

+ 2wJ c Y^k + iA[f73 ' * ]) R ■ (3 ' 19) 
Then by defining the off-diagonal matrix 

Q(x) = ~y f dX to> *( A ' X ) R ( X > *)] ' ( 3 - 2 °) 

J c 

we can compare the integral equation ( |3 .19 ) to the integral equation for the 
matrix Q(x)\&(A, x). As the solution is assumed to be unique we arrive at the 
desired result 

V x +ik{a 3 ,V] =Q^> . (3.21) 
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Property P.8 : We derive now from the above definition of Q{x) the proper- 
ties P. 8 and to that end rewrite ( [3.20 ) as 



q{x) 
r(x) 



p-(X)^ 2 (\x)e 2tXx . 



(3.22) 



C- 



First, for x < 0, we can close the contours with the circle of infinite radius in 
the related half-planes and readily obtain that q and r vanish (remember that 
the contours pass beyond the poles of p). 

To evaluate the boundary values in x — we use ^j^(fc,0) = ^^(A;, 0) = 1 
(the function ^ is continuous in x = 0), but keep the exponentials as they pro- 
duce discontinuous functions in x — 0. Then with help of the large k behaviors 
([T^) and the formulae 



x > 



1 

2i7T 



,~2i\x 



A 



1 

2in 



J2i\x 



A 



= 1 



(3.23) 



we do obtain q(0 + ) = qo and r(0 + ) = ro, where qo and ro are given from the 
spectral data p ± in (^). 

Note that the procedure can be applied to compute the x-derivative of the 
potential by first regu larizing it in x = and by using the next order expansion 
of p(k) given in ( 2.26 ). Indeed we have 



and applying the same procedure we obtain eventually 



-2i\x 



d 



o+ 1 dx 



d 



lim , \ ^~i( x ) r = {v( x ) ~ io&(x)} 



(3.24) 



(3.25) 



The large x vanishing behaviors of q and r follow from their expressions 
( p. 22 ) rewritten as (we consider only q(x), the same computation works as well 
for r(x)) 



q(x) 



1 



+00 -HO 



p + (\)iP+(\,x)e- 2lX * 

n J-oo+iO 



,x e 



-2ik+x 



(3.26) 



71=1 



The above expression vanishes at large x because the path integral vanishes and 
because 



1 



1 



x^oo r+(fc) ' T + (k) 



= , 



(3.27) 



k=k2 



which results from (2.15) and from the very definition of the poles k^. 
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Property P.9 results from the fact that, due to P.l and P.7, the function $ 
also solves Zakharov-Shabat, i.e. 



<Z> x = ik[$,a s }+Q(x)$ . (3.28) 

Then by inserting in it the large k behavior given by P. 5 we must arrive at 
p.9. To check consistency we have then to demonstrate that the following two 
expressions of q(x) 

--/ p+(X)^ 1 (X,x)e-^+2 l Y / Pi^M,x)e-^ , (3.29) 

i /.+oo-i0 N + 

-- / p + mU\ 1 x)e-^ + 2 l Y J P^4M^)z- 2lktx » (3-30) 

do coincide. This simply results from the starting hypothesis of a spectral 
transform p(k) which is continuous and bounded on the real axis (no poles). 
Then both path integral hereabove are equal. 

Note that the relation P.9 giving Q(x) can be evaluated in x = 0, which does 
lead to the expected result thanks to the boundary behavior (2.S) of <i> in x — 
and the large k behaviors of p + and p~ . Note also that, in order to compute 
q(0 + ) one sh ould first deform the contour to pass above all poles k+ and then 
make use of ( 3.23 ). 

3.3 Reductions 

The 3 particular evolutions that we shall study then result from a reduction 
from the two potentials r and q to only one. The first reduction concerns SRS 
and reads 

Q = a 2 Q o- 2 <&r = -q. (3.31) 



Then by direct computation one easily proves that U2^{k, x)G2 solves the same 
spectral equation ( 2.1) as ^>(k, x) and possess the same boundary behavior (2.14) 
in x = 0. Those two functions then coincide, 



*(£;, x) = cr 2 *(fc, x)a 2 , (3.32) 

and examination of their boundary behavior as x — > oo gives the reduction 
constraint on the spectral transform 



T -(k) = T+(k) , p-(k) = -p + (k) . (3.33) 
Considering the poles and residues, the above expressions imply 

N+ = N- , ~p^=-p+, K=k+ . (3.34) 
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The second case is related to SG in thge light cone, it reads 



Q = 02Q01 <=> r = —q 

In that case it is easy to prove that 

$(k,x) = a 2 ^(-k,x)a 2 

and examination of the boundary behavior as x 
reduction constraint on the spectral transform 



P +(M) = -p~(-M) > 

N+ = N~ , p- = -pi 



■(-*) 



and 



■(*) 



(3.35) 

(3.36) 
gives the 

(3.37) 



In the case of sine-Gordon, we also need to implement the reality constraint 
q 6 R for which both relations (3.31) and ( 3.3 5| ) have to hold together. As a 
consequence, by direct manipulations of relations ( 3.33| ) and (3.37), the resulting 
constraint on the spectral transform p = p + which comes in addition to (3.37) 
reads 



p(-k) = p(k) 



(3.38) 



As a consequence the poles in the complex plane come by pair, namely if k n is a 
pole of p(k) with residue p ni then — k n is also a pole with residue — ~p n . In short 



Res{p(fc)} = p n 



Res{p{k)} = -p n 

-kn 



(3.39) 



4 The Dirichlet problem for SRS 

4.1 Boundary values and Lax pair 

The formalism of the direct and inverse problem for the Zakharov-Shabat system 
on the semi-line x > is now applied to the Stimulated Raman Scattering (SRS) 
equations for the pulse envelopes of the pump a(k, x, t), the Stokes b(k, x, t) and 
the medium excitation q{x, t) 

d x a = qbe 2lkx , 8 x b = -qae- 2lkx 

d t q = - J g{k)dk abe- 2lkx , (4.1) 

where g(k) measures the coupling (related to Raman gain) and the parameter 
k is a missmatch wave number resulting from group velocity dispersion Jig] . 
The Lax pair for SRS has been written in the case g(k) = S(k — ko) in ]28| ] 
and regularized by means of the parameter k in [ p0[ which found its physical 
meaning in |l6| ]. 

Dirichlet conditions are prescribed on the line x — , t > 0, 

a{k,x,t)\ x=0 =A(t) , b(k,x,t)\ x=0 = B(t) , (4.2) 
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and on the line x > , t = 0, 



?0M)lt=o = C (X) ■ ( 4 - 3 ) 

Note that we could assume boundary values A(k,t) and B(k,t) instead of A[t) 
and B(t), but here we have chosen to leave the fc-dependence in the arbitrary 
distribution g(k). This is equivalent to assuming separable dependence of A an 
B in k and t. The given boundary value c(x) is assumed to be of Schwartz type 
(it decreases at large x with all its derivatives faster that any polynomial) but 
the physically relevant case is a medium initially at rest, namely c(x) = 0. Note 
that the compatibility of this set of Dirichlet boundary data is infered from the 
system ([O]) which gives 

d x a(k,x,0)\ x=Q =c(0)B(Q) , d x b(k,x,0)\ x=o = -c(0)A(0) . (4.4) 

The above Dirichlet problem has been solved in |ll| and, for sake of com- 
pleteness, we give hereafter the main lines of the derivation of the solution. The 
vector (a,be 2lkx ) T is expanded on the basis i/'i" and (^e 2%kx ^ the coefficients 
being obtained by examination of the values in x = 0. The result reads 



a 



(A - P B)^+ + B0+ . (4.5) 



As a consequence we have 



a \ ( (A- pB)/r 

b J x^oo \ p(A- pB)/r + Bt 



(4.6) 



The time dependence is plug into the spectral data by requiring that the 
solution p + = (V>i j 4>t) °f the Zakharov-Shabat equation be also solution for 
k e M of 

p+ = V+p+ + p+e- ika3X n+e ikaaX , (4.7) 

+ 1_ f g(X)dX f \ a \ 2 - \b\ 2 2abe- 2lXx 
iij X-(k + i0)\ 2abe 2lXx \b\ 2 - \a\ 2 

Note that, as the matrix V(k, x, t) is discontinuous when k crosses the real axis, 
we have to chose a given half-plane. A similar equation could be written for 
p~ = {4>i ,^2) with V~ and Q~ but, due to the reduction symmetry, it is 
redundant. The compatibility between the two operators ( 3.28| )(4.7) with the 



constraint Q x = 0, results in the equation (|4.1|) within the reduction (3.31) 



4.2 Evolution of the spectral transform 

The free entry f2 + = Q + (k,t) is the dispersion relation and, according to the 
general theory ^0|, it is stricly dependent of the choice of the solution. It is 
called dispersion relation because, in the infinite line case, it coincides with 



the dispersion law of the linearized evolution. The differential equation (4.7) is 
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evaluated at x — and x — > oo. This provides 8 equations for the 4 components 
of fi + (fc, t) and the 4 evolutions equations of p(k, t), r(k, t) and their conjugates 
(these evolutions are compatible which actually results from the reduction) . We 
obtain 



n + (*,t) 



2p(k 1 t)ABa 3 



\A\ 2 Z \B\ 2 
2AB 







(4.8) 



(note that the dispersion relation depends on the spectral transform itself), and 
the Riccati evolution 



Pi 



-l[AB p 2 -(\A\ 2 -\B\ 2 ) p-AB] 



(4.9) 



where we have defined 



l(k) = J 



g{X)d\ 
X-k 



Im(fc) > 



(4.10) 



The main fact is that the above evolution conserves the Laurent expansion 
of p(k, t). Indeed, seeking a solution of (4.9) as the power series 



(4.11) 



allows to prove by induction that Pj(t) = for j > 2. For j = 1 we have the 
following evolution equation 



P^= l -AB j dX g{\) = a(t) 



which is solved to get 



Pi(0) 



l-pi(0) f Q a{t) 



(4.12) 



(4.13) 



Then if p\ vanishes at time t = 0, it vanishes for all t. Such is also the case for 
whose evolution (with p\ (t) = 0) simply reads 



Pi 



(0) 







(4.14) 



Consequently, as p\ and p(°> vanish at initial time (because p(k, 0) is the spec- 
tral transform of the Dirichlet datum q(x, 0)), they vanish for all time and the 
solution of the Riccati equation ([0]) possess the right behavior for large k, 
compatible with ( 2.26 ). 

The coefficient of A: -1 gives rise to 



d tP W(t) = ~AB j dXg(X) . 



(4.15) 
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The evolution (4.1 



) written in x = gives 



d t q (t) = -AB J dX g(X) , (4.16) 
and moreover we have, by construction, 

,W(0) = ^. (4.17) 



Consequently the unique solution to (4.15) is 



> (1) W = ^ (4.18) 



as required by (2.26) 



The time evolution (4.9) of the spectral transform has coefficients depending 
only on the input boundary values A(t) and B(t) and contains all the relevant 
information on the nonlinear scattering of the three waves q, a and b. In some 
cases this Riccati equation is explicitely solvable, for instance when A(t) and 
B(t) are proportionnal (which is physically meaningful) as illustrated in 
It can then be checked on the explicit expression of p(k, t) of jl6) that it does 
obey the good large k asymptotics. The fact of being able to provide an explicit 
expression of the spectral transform in a class of physically interesting cases 
makes the very interest of the present approach of boundary value problems. 

4.3 Why it works 

The Dirichlet open-end problem is well posed if the solution q(x,t) constructed 
hereabove is proved to obey for all t the boundary value 

?(*>*) l a =o+ = *>(*) > ( 4 - 19 ) 
and to belong to the class of functions such that 

lim q(x,t) = . (4.20) 

x — >oo 

This proof goes along the following steps. 



Step 1: the solution p(k,t) of the Riccati equation ( J4. 9| ) possess the large k 
expansion 

p(k,t) = ^- + 0(l/k 2 ) (4.21) 



resulting from (4.1q). This result is central to the theory as we see here that 



the time evolution of the spectral transform conserves the properties it has at 



t = 0, sec ( p.26[ ). Note that the expansion can be pursued at any order. 



18 



Step 2: with such a behavior of p(k,t) the property P.3 still holds, namely 



*(M,t)= ( I i) ■ (4-22) 

The proof of the above boundary value is obtained by following exactly the 
method of section [}] where qo has simply to be understood as qo(t). 



Step 3: compute the limit of q(x,t) as x — » + out of the expresion (3.26) by 
following the method of section H used to prove property P.8. 



Step 4: compute the limit of q(x, t) as x — > oo out of the expresion (3.26) 
again by following the proof of property P.8. Note that the essential property 
to be used is the fact that p(k, t) is continuous and bounded for k S R at any 
given t. Particular attention is required when, as t evolves, a pole of p is "seen" 
to cross the real axis. It has been demonstrated in jl8| that the process is 
continuous, or put in other words, a resonance (pole in the lower half-plane) 
dissapears and a pole appears while never belonging to the real axis. The same 
situation arises for the Schrodinger eigenvalue problem |p5| . 

4.4 Summary of the method of solution of SRS 

The solution of SRS with Dirichlet condition is achieved along the following 
steps: 

1 - The spectral transform at time zero p(k, 0) is obtained form the Dirichlet 
condition q(x) by finding the vector solution ^{k, x, 0) of the integral equations 



(2.4) in the upper half fc-plane, namely 

POO 

0+(fc,x,O) - / dy q(y)<p+2(k,y,0)e- 2ik ^ 



0+ (k,x,0) = 1 - f dy q(y)<fa(k,v,Q) , (4.23) 
Jo 

and then p(k,0) results from 

p(M) = # 2 (M,0)| a=0+ • (4-24) 

2 - p(k, t) is calculated in the upper half /c-plane from p(k, 0) by solving the 
Riccati time-evolution 

p t = ~l[AB p 2 -(\A\ 2 -\B\ 2 ) p-AB] , (4.25) 

with entries A(t) and B(t). The N poles k n in the upper half-plane and related 
residues p n of p(k, t) must be identified for each value of the parameter t. 
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3 - The integral system (3.5)( |3.6| ) which, within the reduction ( |3.33 ), sim- 
plifies to 

f/>ti(Js,x,t) =1 + 77- T^r p(^t)e 2iXx ^ 2 (\x,t) 

N(t) 



J2 =^^J^{K{t),x,t)e 2i ^ , (4.26) 



-i p+oo+%0 j\ 

=— / — — p(\,t)e- 2i ^+(\x,t) 

2«7T J-oo+iO X ~ k 
N(t) . . 

E klih^(k n ,x,t)e- 2i ^ , (4.27) 

71=1 

is solved to get ip^^k, x,t). This system of equation holds for any fixed value 
of the "parameters" x and t. 



4 - Finally the solution q(x,t) is obtained from expression (3.26), i.e. 

r + oo+iQ 



q(x, *)=--/ dX p(X, t)il>ti ( A > t)e- 21 v 



N(t) 



+ 2i^2p n (t)^+(k n (t),x,t)e- 2ik ^ x . (4.28) 

n=l 

5 The Dirichlet problem for sine-Gordon 

5.1 Boundary values and Lax pair 

Another case where the method applies sucessfully is the sine-Gordon equation 
in the light cone 

xt + sin0 = O, (5.1) 

which is solved here for Dirichlet conditions 6o(t) and 9{x), namely 

te[0,T] : 0(0, i) = 9 {t) , xe[0,oo) : 9(x,0)=6(x), (5.2) 

in the class of functions vanishing (modulo 2ir) with all derivatives as x — > oo 
for every value of time t (the value T is arbitrary) . 

This equation is fundamental in many areas of physics and has been con- 
sequently widely studied. The solution of the boundary value problem on the 
quadrant also has been the subject of intense studies, see e.g. J|, ||, 10, [b| [u| 
and the references therein. 
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The Lax pair, found in J29|, is given by the Zakharov-Shabat equation for a 
particular choice of potentials r and q, namely 

$ x = ifc[$,(T 3 ]+Q$ , Q = ~a a e x , (5.3) 
together with the following evolution 

$f = - 4r (<r 3 cos 6 + (7i sin 0) $ + $ e - lfc<T3a W fc,T3a: . (5.4) 
4fc 

We have chosen to write this Lax pair for the particular solution $ defined in 
section |[ Note that t he ab ove particular structure of the potential Q obeys the 
reduction constraint ( |3 .35 ). 



5.2 Evolution of the spectral transform 

By using the behaviors fl2.8|) and (|2.9|) in (|5.4|) we calculate f2 and get the time 



evolution of the spectral transform. Within the reduction relations (|3.37| ) we 
eventually obtain 



cos 6*o — p(— k,t) sm0 o \ , > 

ifc^ - cos 9 + P (k,t) sin 9 ) ' 1 j 



pt = — [p 2 sin O — 2pcos# — sin# l . (5.6) 
ik 

Note that the reduction constraint ~p(—k) = p(k) which ensures real value po- 
tential 9(x,t) is conserved by the above time evolution (as soon as 9q £ K of 
course). 

The evolution of p(k,t) is of Riccati type and depends exclusively on the 
boundary value 9o(t), while the other Dirichlet condition 0(x) fixes the initial 
value p(k,0). Following exactly the same method as for SRS we show that, as 
Pi(0) and p(°)(0) vanish (by construction), the large k asymptotic behavior of 
the solution p{k, t) reads 

p(k,t) = jp^ +0(l/k 2 ) , d t pM = ~sm0 Q . (5.7) 
k 4 

Thanks to the sine-Gordon equation evaluated in x — we can integrate the 
above time evolution and obtain 

p«( t )=-i W ) = ^, (5.8) 
as indeed, by construction of p(fc,0) from 9(x), we have 



P (1) (0) = --#* = -~0 x (O,O). (5.9) 



1 K 



1 
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The last point to consider is the apparent singularity of p(k,t) when k — > 
by real values. Actually for k — > the spectral transform possess the following 
regular behavior 

p(k, t) = r (t) + kn(t) + 0(k 2 ) , (5.10) 

that can be deduced from the very definition ( 2.1 0| ) of p(fc, t) (or equivalently by 
( |5.18 ) below), by expanding these expressions in power of k. Inserting the above 
expansion in ( |5.6| ) gives the following unique solution vanishing for 9q = 0, 

p(0,t) = r (t) = ^fl. (5.11) 
sin fro 

The other terms r n (t) in the expansion are then recursively given in terms of 
Tj, (J < n), without constraint. 

5.3 Why it works 

The Dirichlet open-end problem is well posed if the solution 9(x,t) constructed 
hereabove is proved to obey for all t the boundary value 

6(x,t)\ x=0+ =9 (t) , (5.12) 

and to belong to the class of functions such that 

lim 9Jx,t) = . (5.13) 

x — >oo 

This proof goes along the following steps. 



Step 1: the solution p(k,t) of the Riccati equation (|5.6| ) possess the large k 
expansion 

P(M) = ^ + 0(1A2) (5-14) 



where we have 



1 , . 1 '* 



qo(t) = ~2^(°) + 2 / s[n9 ° ^ 5 - 15 ) 



resulting from (5.7). 



Step 2: then property P.3 holds for all t, namely 

*(M,i) = ( J J ) ■ ( 5 - 16 ) 

The proof is obtained by following the method of section [| where go has simply 
to be understood as qo(t), given from Dirichlet conditions by ( 5.15| ). 
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Step 3: the limit ofq(x, t) as x — » + is computed out of the expression (3.26) 
by following the method of section || used to prove property P. 8. 

Step 4: the limit of q = —6 x /2 as x — > oo is computed out of the expression 
( p. 26 ) again by following the proof of property P.8. Note that the essential 
property to be used is the fact that p(k, t) is continuous and bounded for k G R 
at any given t. 

5.4 Summary of the method of solution of sine-Gordon 

The solution of sine-Gordon with Dirichlet conditions is achieved along the 
following steps: 

1 - The spectral transform at time zero, p(k, 0), is obtained form the Dirichlet 
condition 0(x) by finding the vector solution </> 2 (k, x, 0) of the integral equations 
wA) in the upper half fc-plane, namely by solving 



cf>Uk,x,0) = dy q(y)^2(k,y,0)e- 2lk(x - y) 

J X 

0+(fc,x,O) = l+ / dy q{y)4>t 2 {k,y,G) , (5.17) 



where q = —9 x /2. Then p(k, 0) results from 

p(k,0)= 4>t 2 (k,x,0)\ x=0+ . (5.18) 

2 - The spectral transform at time t, p(k, t), is calculated in the upper half k- 
planc from p(k,0) and Dirichlet condition 6o(t) by solving for k ^ the Riccati 
time-evolution 

Pt = TT \p 2 sin - 2p cos 9 - sin 9 ] , (5.19) 

and for k = 

K0,t) = ^fl. (5.20) 
sm 0o 

The N poles k n in the upper half-plane and related residues p n of p(k,t) must 
be identified for each value of the parameter t. 

3 - The integral system (3.5)( |3.6| ) which, within the reduction ( |3.37 ), sim- 
plifies to 

Tjj+(k,x,t) =1 - — / — — p{-\t)e 2 ^^ 2 (\x,t) 

N(t) 



E ^fe(^n(t),M)^ Mt)l , (5.21) 

n— 1 n ^ ' 
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1>n(k,x,t) =— / — -r p(A,t)e- 2lA ^+(A,x,t) 

N(t) . . 

- E k^k^ {k ^ x ^ e ~ 2lkn(t)x ' (5 - 22) 

is solved to get tp^^k, X,t). This system of equation holds for any fixed value 
of the "parameters" x and t. 

4 - Finally the solution 9(x, t) is obtained under the form 6 x (x, t) = —2q(x, t) 
where q(x,t) results from tp^ and p as 



1 



7T 



+ oo+iO 



g(x, t) = / dX p(X,t)ipi 1 (X,x,t)e 



-2i\x 



oo+iO 
N(t) 

+ 2iJ2pnmtMt),x,t)e- 2ik «M x . (5.23) 

n=l 

5.5 An explicit example 

To illustrate the above method we consider the explicitely solvable case of a 
piecewise constant boundary datum 9o(t), namely 



t e [0,*i] 
t £ [tj,t j+ i] 

t e [t n , OO] 



%{t) = o , 

%{t) = Vj , (5.24) 
%(t) = , 



compatible with a vanishing initial condition 

6(x,0) = 6(x) = . (5.25) 

Note that such a boundary datum represents an arbitrary localized function 
&o(t) discretized with n-step time grid. 
We adopt the following notations 

pj{t) = p{k.t) for t G [tj,t j+1 ] , 

Pi = Pi(tj+i) = P( k i tj+i) • ( 5 - 26 ) 



Then the solution of the evolution ( |5.6| ) can be written for t £ [tj,tj+i] 
n M _ b j (p j - 1 +a j )e i W* k -a j (p j - 1 +b j ) 

with the following definitions 

1 + b IzSEm. (5.28) 

sin ifj sin ipj 



24 



Note that the limit tpj — *■ out of ( |5.27 ) does produces the solution to evolution 
( |5~G| ) where 60 — ipj is set to zero. 

The large-k behaviour of solution (5.27) is obtain by induction and reads 

1 , 



Pj(t) : 

Hence we have 



1 

4k 



3 

(* - tj) sin <pj +J2(te 



s-ilsin ipe-i 



0( 



k 2 ' 



i 1 



(5.29) 



(5.30) 



compatible with (5.7) and solution (5.27) has the expected large-k behaviour. 
At the last step, that is for (p n = 0, (5.27) can be writen 



Pn(t) = - 



(&„p„_i - l) e *(*-*n)/2fc _ a nPn - X + 1 



(p„_x + a „) e i(t-tn)/2fe _ ( pn _ 1 + bn ) 
which eventually simplifies to 



(5.31) 



(5.32) 



As a consequence, the spectral transform evaluated after the total effect of the 
boundary datum 6o(t), i.e. for t > t n , is completely determined by p n -\ — 
p(k,t n ). This function of k is itself obtained as the following mapping 



Po = , 



Pi 



(pj-i + a 7 -)e i fe+i-*j)/ 2 * ! - ( Pj _! + bj) 



,3 = 1 



1 . (5.33) 



This is an explicit expression that allows to get insight in the effect of a generic 
input datum 9o(t) on a vanishing background under suitable (arbitrary) dis- 
cretization. We report the study of the above mapping (in the complex fc-planc) 
to forthcoming work and simply mention here that the main point concerns the 
study of the creation (anihilation) of poles of p(k,t) in the upper half fc-planc. 

Remembering the reality constraint (3.3£), poles of p(k) come by pair, and 
a pole k — k n (t) of pj(t) is given by 



e 2k n = 



Pj-i + h i 

pj-l + Oj 



(5.34) 



where Pj-i is itself a fuction of k n . 

This equation for the pole k n is in general highly complicated but if Pj-\ 
(as for the initial stage t = t\), the above relation becomes 



2k n — 



1 



COS if 1 



1 + COS ipi 



(5.35) 
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Defining 



In 



the solution of (5.3S) can be written 



1 — COS ifi 
1 + COS (fx 



*»(«) = ^[(2»+i> + *r] (2 J^ + 72 



(5.36) 



(5.37) 



where n assumes values in Z. The above expression is identical to the one 
obtained in Jl^] (see formula 5.52) for t — t\=T, A property of our approach is 
to explictely display the time dynamics of the spectral transform (such as poles 
generation and motion). 





0.5 n 






0.4 






0.3- 






0.2- 






0.1 












Figure 1: Positions of the poles of p(k,t) for n 6 [—20,19], t — t\ = 10 and 
ipi = 0.657T 

Among the above set of poles, only those with positive imaginary part con- 
tribute to the soliton spectrum. It is then remarkable that 

<Pi e [0, \[ =>• 7 < =>• Im(fc„) < , (5.38) 
fx e]|,7r] =► 7 > =► Im(fc„) > , (5.39) 

which shows that no poles are generated by a consant boundary datum 8q (t) — 
ifi smaller than 7r/2, while an infinity of them are generated for ip± > 7r/2. An 
example of poles locations is shown on the figure where we see that the point 
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k = is an accumulation point for the infinity of poles. Indeed when n — > oo 



we have from (5.37) k n — ► 0. The case ip± — tt/2 is singular as p(k,t) would 
then possess poles on the real axis, for which the theory fails. 

Note that the discretization of a smooth localized boundary datum &o(t) 
produces small yjj's. Hence the soliton generation would then result from the 



contributions in equation (5.34) of the Pj-e- In the case of a single constant 
boundary datun of amplitude greter than 7r/2, the assignation of solitons to the 
presence of an infinity of poles is not clear (usually an infinity of poles is related 
to a self-similar structure). 
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